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Abstract 

We employ the field theoretic approach to study the quantum noise problem in the mirror-field 
system, where a perfectly reflecting mirror is illuminated by a single-mode coherent state of the 
massless scalar field. The associated radiation pressure is described by a surface integral of the 
stress-tensor of the field. The read-out field is measured by a monople detector, form which the 
effective distance between the detector and mirror can be obtained. In the slow-motion limit of 
the mirror, we are able to identify various sources of quantum noise that lead to uncertainty of the 
read-out measurement. Since the mirror is driven by radiation pressure, the sources of noise, other 
than the shot nose given by the intrinsic fluctuations of the incident state, may also result from 
random motion of mirror due to radiation pressure fluctuations and from modified field fluctuations 
induced by the displacement of the mirror. Correlation between different sources of noise can be 
established as the consequence of interference between the incident field and the reflected field out 
of the mirror in the read-out measurement. The overall uncertainty is found to decrease (increase) 
due to the negative (positive) correlation. In the case of negative correlation, the uncertainty can 
be lowered than the value predicted by the standard quantum limit. We also examine the validity 
of the particle number approach, which is often used in quantum optics, and compared its results 
with those given by the field theoretical approach. Finally we discuss the backreaction effects, 
induced by the radiation pressure, that alter the dynamics of the mean displacement of the mirror, 
and we argue this backreaction can be ignored for a slowly moving mirror. 
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I. INTRODUCTION 



The measurement of a very weak force such as the detection of gravitational wave demands 
a very high resolution interferometer. The role of the interferometer is to transform small 
variation of the relative spatial separation of the mirrors between two arms into variation 
of the output photons. Quantum noise presumably places a limitation on accuracy of the 
measurement [1]. 

In quantum optics, the discussion related to quantum noise of radiation fields on the 
mirror has been mostly based on the photon-number approach. It is argued by Caves [3l H] 
that quantum noise in a laser interferometer may come from the quantum nature of the 
light directly via photon number fluctuations (shot noise) or indirectly via random motion 
of the mirror due to a fluctuating force (radiation pressure fluctuations). Minimizing the 
total uncertainty from two sources of noise, with an assumption of no correlation between 
them, may give the standard quantum limit (SQL), when an input power of the light is 
appropriately tuned. Additionally, Caves further suggests that injecting squeezed vacuum 
state into the unused port of an interferometer should lead to the SQL with a lower optimum 
input power [3], as compared with the standard interferometer ||4j. This seminal work also 
shows the manipulated quantum fluctuations of the electromagnetic fields are responsible 
for both the radiation-pressure forces and the fluctuations of the photon number. Then, 
an inquiry into the assumption of no correlation between these two noise sources arises. In 
particular, the presence of negative correlations may lead to the uncertainly even below the 
SQL. To establish their correlations requires an unified scheme that has been proposed by 
Loudon [7J. The judicious use of the correlated squeezed state is proposed to establish the 
above correlations with which to possibly push the sensitivity of the interferometer beyond 
the SQL [inj. The redesign of an interferometer by adding some optical components to either 
manipulate the read-outs or modify the dynamics of the test mass, such as in the "optical 
bar" [TT] or the "optical spring" [12] scheme, is also possibly to give non-zero correlations 
in order to beat the SQL. 

The coupling between the mirror and the radiation field can be obtained by considering 
the mirror as a reflector that imposes boundary conditions of the field on the surface of the 
mirror. In 3+1 dimensions, field quantization is treated perturbatively for a slowly moving, 
perfectly reflecting mirror by Vilenkin and Ford [T3], with the boundary conditions of the 
field on the surface of the mirror at all times. They then explore the backreaction dissipative 
effect on a moving mirror, induced by vacuum fluctuations of the field. Later a corresponding 
Langevin equation was derived by two of us [H] with the method of influence functional, in 
which the accompanying noise, manifested from vacuum fluctuations is taken into account. 
Thus, the effects of the backreaction dissipation and its accompanying noise can be related by 
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the underlying fluctuation-dissipation relation. Unruh [8] extended the quantization scheme 
to the interferometer case where the field propagates toward a mirror and then acquires a 
time- dependent phase shift after reflecting off the mirror. Because of the moving boundary, 
the phase shift in general depends on the mirror's motion, but can be obtained by assuming 
slow-motion of the mirror. Similar phase shift was also introduced by Kimble et al. [9], 
and it has been proposed to improve performance of the faster interferometer, particularly, 
in [ig. 

Here we adopt the field quantization scheme similar to Unruh's proposal to investigate the 
problem of quantum noise and its associated correlation in a simpler configuration, shown 
in Fig. [TJ The basic idea is to consider a free, perfectly reflecting mirror. A coherent state 
of the quantum scalar field is normally shined on the mirror's surface, and exerts a pressure 
force on the mirror. It drives the mirror into motion, which in turn leads to variation of 
the reflected field. Then the read-out field is measured by a standard monopole detector, 
which is placed somewhere between the mirror and the field source. Since the radiation 
pressure is described by a surface integral of the stress-tensor of the scalar field, if the input 
quantum state is not an eigenstate of the operator associated with the response function of 
the detector or the stress-tensor, the resulting measurement will exhibit fluctuations. 

The goal of this paper is two-fold. First of all, we lay out a field-theoretic approach by 
which the effective distance between the detector and the mirror can be decoded from the 
read-out fields. In this simpler system, the sources of the uncertainty in the read-outs can 
be easily identified in the slow motion limit when the mirror is driven by radiation pressure 
of the incident field. Thus reduction in uncertainty of this effective distance is of essence to 
improve the sensitivity of the interferometer. When the incident field is in the single-mode 
coherent state, it will be shown that the particle number approach, typically used in quantum 
optics, is equivalent to this field theoretic approach in the late-time limit with a large number 
of the scalar particles (photon). We will elaborate the previously mentioned approximations 
later. In addition to shot noise and radiation pressure fluctuations, a new source of noise 
comes from modified field fluctuations. It is induced by the mean displacement of the mirror 
as a result of the radiation pressure. Its effect can be systematically examined in the field 
theoretic approach. Secondly, the emphasis is put on how correlation among all sources 
of quantum noise associated with the radiation field, can emerge when the interference of 
the incoming field and the reflected filed by the mirror are taken into account. We further 
examine the effects of correlation on possible reduction of uncertainty in the read-outs of the 
monopole detector. We show an example in which the overall uncertainty can be suppressed 
to beat the standard quantum limit. 

Our presentation is organized as follows. In Sec. |TT], we introduce the quantization scheme 
of a massless scalar field subject to boundary conditions on a free, perfectly reflecting mirror 
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when it undergoes slow motion. The motion of the mirror arises from radiation pressure of 
the scalar field. A monopole detector is introduced to measure the readout field. We then 



identify the sources of quantum noise from the fluctuations of the readouts in Sec. Ill with 
an emphasis on establishing correction between them. The effects of correlation are further 



studied in Sec. |lV]for possibility of reducing overall quantum noise. In Sec. |V]the Langevin 
equation of the mirror, including backreaction from the radiation field, is derived, and we 



show that this effect is negligible for slow motion. Conclusions are drawn in Sec VI 



The Lorentz-Heaviside units and the convention h = c = 1 will be used unless otherwise 
mentioned. The signature of the metric is diag{?7^y} = (+1, —1, —1, —1). 



II. FIELD QUANTIZATION AND DETECTOR THEORY 

We consider a free mirror of perfect reflectance is illuminated by a single-mode coherent 
laser, which propagates in a direction normal to the mirror. The radiation pressure of 
the laser will drive the mirror into motion, and the field, bounced off the mirror, can be 
measured by a detector. For simplicity and without loss of generalization, the radiation 
field is represented by a massless scalar field (p to suppress complexity caused by the vector 
nature of the electromagnetic field. The mirror, with mass m, is placed at z = L. Its lateral 
dimensions in the x, y directions are large compared to wavelength of the incoming state, 
and the cross sectional area is A\\. We use xy to denote directions normal to the z axis. The 
incident field is a plane wave propagating along the positive z direction toward the mirror. 
The boundary condition of the field on the mirror surface is then given by 

0(x||,L + g(t),t) =0, 

where q{t) is the displacement of the mirror from its original position ai z = L due to 
radiation pressure. The approximate solution to the field equation subject to this moving 
boundary condition in the region of 2; < L is [S] 

= 0+ + 0- , 

where the positive-frequency solution 0+ is given by 

/' d^k 1 
(^^^«fc^k(x,i), (1) 

with 

The negative-frequency solution 0~ is the Hermitian conjugate of the positive-frequency 
solution 0"*", that is, 0~ = (0"^)^- The prime over the integrals of k means that the modes 
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with kz < are excluded because these modes propagate in the negative z direction and 
will not interact with the mirror in the z < L region. Here we have assumed that the mirror 
undergoes slow motion (g ^ 1). Since the incoming wave is a right-moving plane wave along 
the z direction, the left-moving reflected wave will receive a time-dependent phase shift due 
to the motion of the mirror. The phase shift of the reflected field at time t and position z 
depends on the mirror's position at an earlier time Ir = t ~ \L + q{tR) — z\ if we take finite 
lapse of propagation into consideration. When the displacement is small, we may take the 
retarded time a,s = t ~ \L ~ z\ in the leading order approximation. 

Due to small displacement approximation, the force acting on the mirror by the fields can 
be approximated by the area integral of the zz component of the energy-momentum stress 
tensor of the scalar field, evaluated at the mirror's original position z = L, 

Fo(t) = y"d^||T,,(x||,L,t), (3) 

where the zz component of the energy-momentum stress tensor takes the form 

TU^\l,z,t) = \ (9,0)'- . (4) 

Hereafter we will use dA\\ to denote the area element normal to the z direction. Thus, if 
the mirror is at rest aX z = L for t < 0, then the displacement operator due to the action of 
radiation pressure of the incident state is formally given by 

q{t) = - fds fds' [ dA\\ T,,(x||,L,s') . (5) 
Jo Jo J 

Since the energy-momentum stress tensor is quadratic in the field operator 0, the displace- 
ment operator is not well-defined. A proper renormalization of mirror's position operator 
can be done by absorbing the vacuum expectation value of the displacement operator into 
the initial position L, namely, 

q{t) + L = q{t) - (g(t))o + {L + (g(t))o) = g.(t) + U . 

This prescription is physically sensible since when the mirror is placed at its initial position, 
the vacuum fluctuations of the field have already interacted with the mirror before the 
incoming state impinges on it. Thus, the renormalized displacement operator qr{t) now 
takes the form 



qr{t) = — ds ds' dA\\ |T,,(xj|, L, s') - (T,,(x||, L, s'))c 

Jo Jo J 

Here the backreactions of the field is ignored. This backreaction can be incorporated in 



the Langevin equation, Eq. (43), and is found negligible for a slowly moving mirror. To 
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FIG. 1: Schematic diagram of the field-mirror system. 



avoid cluttering the notion, we will suppress the subscript r of the renormalized quantities 
hereafter. 

To measure the outcomes of the scalar field, we employ a monopole detector. We assume 
that the detection processes is based on the stimulated transition of the detector due to 
coupling between the scalar field and the monopole moment m{t) of the detector, 

"t 

ds m(s) 0(x, s) . (6) 



In the standard detector theory, the dynamics of the detector is governed by some time- 
independent Hamiltonian Hd. Assume that the detector is a two- level system with eigen- 
energies Ei, E2 such that He,\Ei^2) = -E'i,2|-£'i,2) and Ei < E2. In the interaction picture the 
monopole moment evolves in time as 



m(t) 



m(0) e 



(7) 



Then by first-order perturbation theory, the transition rate between these two states of the 
detector is given by 



P(Ei ^ E2) = |(^i|m(t)|E2)P X Ii^{E2 - El) . 
The response function II,^ is defined by 

n^(E) = [ds [ ds' e"^^'-''^ (0-(x, s)0+(x, s'))t 



(8) 



5{E — u) j ds {(f) (x, s)0^(x, 



(9) 



where we have assumed that the incident field is in a single- mode coherent state \a) with 
frequency u. Eq. (|9| implies that we may construct an integral quantity by 



I{zo,t; L + q) = J ds j dA\\ S{-x\\, zq, s; L + q) 
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with 

S{x; L + q) — (f)~{x) (l)'^{x) , and x — {t,x\\, zq) , (10) 

where the detector is placed at z — zq and the position of the mirror is located at 2; = L+q{t). 
The quantity / is a measure that generically summarizes accumulated measurement results 
over finite time duration. Integration over an area of the mirror size illustrates the situation 
that if the size of the detector is much greater than 27r/k\\ of the incoming coherent mode, 
then the measurements have been averaged in the lateral spatial direction. Then in terms 
of the mode functions of the scalar field, I{zo, t;L + q) can be expressed as 

I{zQ,t;L + q)= I ds I dA\\ j -—^ I -—^ ^=^=Ul{^\\^ 
Jo J J (27r) 2 J (27r) 2 \/2uj \/2u}' 

(11) 

For the coupled mirror-field system, we notice that they are at least two sources that will 
contribute to the uncertainty associated with I. One directly results from quantum fiuctu- 
ations of the incoming field, but the other indirectly comes from the stochastic motion of 
the mirror, which in turn is the consequence of radiation-pressure fiuctuations of the field. 
In other words, the quantum fiuctuations of the incoming field will drive the mirror into 
stochastic motion, and next the moving mirror will modify the refiected field in a compli- 
cated way, depending on mirror's motion. If we use the a detector to measure field variation 
at 2; = ^0 < the detector will see combined effects from the incoming and the refiected 
field. It will be shown later that the interference between the right-moving (incoming) and 
left-moving (outgoing) field leads to the result that the phase of the expectation value of the 
operator I{zq, t; L) depends on the effective distance between the mirror and the detector. In 
addition, this interference effect is essential to establish correlation between different sources 
of quantum noise. This phenomenon gives a hope to reduce overall quantum noise. 



III. QUANTUM NOISE 

Since the mirror is driven by the radiation pressure of a single-mode coherent state with 
the frequency a), the assumption of slow motion of the mirror can be realized by 

uq{t)^l, q{t)/L<^l. (12) 

This suggests that wc may formally expand the field operators, 0^(a;; q+L) and I{zq, t; q+L), 
in a power series of the the displacement operator q{t). The operator ordering can be quite 
ambiguous in this case. However, it can be shown that the leading-order result is independent 
of operator ordering. Thus, for example, up to the second order in q, we may just expand / 
into 

I{zo, t;q + L)^ I{zo, t; L) + q{t,,,) dJizo, ^'^^) + \ ^'(^-t) ^^/(^o, t; L) + e»(?=^) , (13) 
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irrespective of the operator ordering. Apparently, the first term represents the measurement 
performed with respect to the mirror's original position. The other terms can be understood 
as corrections which result from the mirror's motion due to radiation pressure. 

The single-mode coherent state of the background field can be obtained by applying the 
displacement operator D{a) on the corresponding vacuum state, 

\a) = D{a)\0), (14) 

where the displacement operator D{d) is defined by 

D{a) = exp [a — a* a^] . 

Then, for a fixed position of the mirror, the signal is given by the expectation value of 



' I{zo,t; L) )a, which is obtained from Eq. (11) by setting q = 



{I{zo,t;L) )a = 77— - sin^w(L - zq) x t . (15) 

The result depends on the relative distance between the detector plane and the mirror 
surface, which arises from the optical path difference between the right-moving (incoming) 
field and the left-moving (outgoing) field. If a free mirror, instead of a fixed one, is considered, 
the above result will be modified to account for the additional effect due to the mean 
displacement of the mirror. We further assume that the coherent state under consideration 
has a sufficiently large particle density lap so that in this state, the ratio of the operator's 
fluctuations over its expectation value, which is typically inversely proportional to lap, is 
small. Thus, the variation of I{zo, t;q + L) can be approximated by 

AI{zo, t-q + L) = AI{zo, t; L) + { ddizo, t; L) )^ Aq{t,et) + ( g(tret) )« AdJizo, t; L) 

+ ( g(tret) )a ( dll{zo, t; L) )^ Aq{t,,,) + ^ ( g(t,et) )l Adll{zo, t; L) 

+ O(A') , (16) 

if, given an operator O, its variation AO is defined as AO = O — {0)a- This gives a 
leading-order effect of the noise. Later we will offer an order-of-magnitude estimate of the 



ignored terms. The first term of (16) arises from the intrinsic fluctuations of the incident 
state. For a large value of |ap, the particle number has relatively small uncertainty; thus 



according to (15), the main contribution to the variation of I{zQ,t;L) for fixed uj does not 



come from amplitude uncertainty, but is given by the phase fluctuations, as will be seen 



later in Eq. (30). It is thus associated with shot noise. The other terms may be understood 
as the induced effects due to either the position fluctuations of the mirror, or the modified 
field fluctuations as a result of the mean displacement of the mirror. Then the variation of 
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I{zq, t;q + L) allows us to define the overall uncertainty of the effective distance between the 
detector and the mirror, 

^ AI{zo,t-q + L) 
\{dLl{zo,t;L))^\' 

as long as | uAz \ < 1. The normalization factor ( OlIIzq, t; L) )g is to measure the change of 
1(2:0, t) due to variation of the mirror's position. Control of the above uncertainty is essential 
to improve the sensitivity of the interferometer, in which the measurement of the separation 
between the mirrors and the beam splitter plays an important role. 

To compute the square of the effective displacement uncertainty (A^;)^ of the mir- 
ror, we may encounter terms such as {S{x) S{x')) — {S{x)) {S{x')), {Tzz{x) T^zix') ) — 
{Tzz{x)) {Tzz{x')), and their cross correlation terms. Since both S{x) and Tzz{x) are 
quadratic in field operators they are not well-defined in the coincidence limit x — x'. 
We expect that the square of the effective displacement uncertainty should contain terms in- 
volving a product of four field operators. Thus it may be more transparent if we decompose 
such a product as follows, 

01020304 =: 01020304 : + : 0102 : (0304)o+ ■ 0103 : ( 0204 )o 

+ : 0104 : ( 0203 )o + : 0203 : ( 0104 )o + : 0204 : ( 0103 )o 

+ : 0304 : ( 0102 )o + ( 0102 )o( 0304 )o + ( 0103 )o( 0204 )o 

+ (0104 )o( 0203 )0, (18) 

where (■ ■ ■ )o denotes the expectation value of the scalar field in the Minkowski vacuum. The 
first term is a fully normal-ordered term, the next six terms are cross terms, and the final 
three terms are pure vacuum terms. For a single-mode coherent state, the fully normal- 
ordered terms cancel, and then we are left with 

( 01020304 )a " ( 0102 )a ( 0304 )a = {■ 0103 : )a { 0204 )o + Othcr CrOSS tcrmS 

+ ( 0102 )o ( 0304 )o + other pure vacuum terms . 

The cross terms and the pure vacuum terms may be singular because the field operators 
will be evaluated at the same point. It has been shown that the spacetime average of 
the pure vacuum term leads to a result that varies as an inverse power of the size of the 
spacetime [I7j. On contrary, the cross terms depend on the particle density of the coherent 
state. For a large value of the particle density, the cross terms can dominate over the pure 
vacuum terms [T7] . Thus we will consider the contributions from the cross terms and ignore 
the pure vacuum terms. Since the integrals that contain cross terms may still diverge, a 
finite result is obtained [17] by choosing the principal value of the singular integral. Here 
we take an alternative approach to cope with the singular vacuum contribution so that all 
results become regular. They prove to be the same in the end. 
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It is straightforward to compute the expectation value of the operators in Eq. (16). Let 
us focus on the large t ^ L ~ zq limit. They are given by 

^11 



\a\' 



X 



(19) 



(27r)3m 

The results of {dLl)a and can be obtained by taking the derivative of the expres- 

sion Eq. (15) with respect to L. 



Let the dispersion associated with Az be defined as 

{^P{zo,t-q + L))o, 



{Az') = {{AzY) 



[dLlT{zo,t;L)) 



(20) 



From (16), the whole expressions in (A^;^) can be grouped according to their physical mean- 
ing. They respectively come from the shot noise (sn) contribution associated with intrinsic 
fluctuations of the incident coherent fields, and the contributios of radiation pressure fluctu- 
ations (rp) and modified field fluctuations (mf), both of which are induced by the mirror's 
motion, 



{Az'^)mf 



{dd{zo,t-L))l ' 
(Ag^(t))., 



2 {A{dLl)\zo,t-L))^ 



{dLl{zo,t;L)) 



(21) 
(22) 
(23) 



In addition and more importantly, there exist the cross terms owing to correlation between 
different sources of uncertainty. Up to the order g^, we lump these terms together in 



1 



{{AI{t),Aq{t)})^ + 



'dLmh 



{AdLl{t),Aq{t)})^ + 



+ 



'At)) 



2{dLl{t))l 



{{Adll{t),AI{t)})^. 



{{AI{t),AdLl{t)})^ 

{A/(t),Ag(t)})^ 
(24) 



{dLl{t))l 

q{t))^{dim) 



From now afterwards, the quantity I and its derivatives are understood to be evaluated at 
z = L, and we replace I{zo, t; L) simply by its shorthand notation I(t). We will discuss these 
contributions individually in the following sections. In particular, the emphasis will be put 
on the correlation effects that may give the hope to reduce the overall uncertainty. Here we 



note that Eqs (21), (22), (23), (24) involve the expectation values of the anticommutators 
of the operators, so in the leading-order approximation we will obtain the same results 



irrespective of operator ordering on expanding Eq. (11). 
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A. shot-noise terms 



To the overall dispersion of A2;, the contribution from the intrinsic fluctuations of the 
incident field is given by the term ( AJ^ )a, which is expressed as 



ds / dA 



(27r)ti (27r)fyo 



H^o,t; L))l 

' d^k r d^k' 



ds / c/^jj 



(27r)2i (27r)i^ ^2k ^2k' ^2k V2k' 



X Wk(x|j , zo, s; L) Wk' (x|| , zq, s; L) W?(xj| , zq, , s'; L) Wj^, (xy , zq, s'; L) 



X 



^1 



(2vr) 



q;| — sin — zq)] j ds I ds 



duj 1 
27r2w 



where the area integral is carried out over a measuring plane at z = zq, and its total area A\\ 
is assume to be large but finite. Performing the time integrations and using the fact that 



we have 



lim — 

t-s>oo CO — OJ 



Al\zo,t;L)), 



sm 



(CJ — ljj)t 



All 



i2n) 



\a\ 



'k6{u — COq) 



sin''[w(L — zq)] X t . 



(25) 



(26) 



Together with the normalization factor {8^1)^, Eq. (21) gives the positive contribution 

1 1 



tan [ijj{L — zq)] 



{dLl)l Put A 
where the energy flux P of the incident state is defined by 

Ai 



t>l/a;, and t > L - 2:0 , (27) 



P 



i2ny 



(28) 



The result (27) is inversely proportional to the power P, and that is the typical behavior of 
the shot nose. 

Since the read-out field I{zo, t; L) is to count the number of the particles being detected, 
the same result can be obtained in terms of the particle number operator n if we make such 
an identification, 



/(zo, t; L) = — sin [ui{L — zq)] x n 



UJ 



(29) 



and notice that for a coherent state, {n'^)a — {n)"^ = {n)a, where the mean number of 
particles {n)^ that reaches the surface A\\ within time t is that {n)a 



^11 

(27r)3 



a H. This 



shot noise can be understood as phase fluctuations of the field [26]. The phase uncertainty 
of I can be effectively defined as 



A0 



(30) 
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The variance of the phase A0 is found proportional to 1/P as well. It is worth to mention 
that this result depends on not only the input power P but also the distance between the 
mirror and detector due to interference between the incident and reflected fields. The proper 
choice of the parameter L — zq can be used to avoid the potentially large uncertainty due to 
the factor tan^[a;(L — zq)] and then to minimize the overall uncertainty. The same strategy 
has been used by Caves in [3l H] to finding the appropriate spot in the fringe pattern so as 
to further reduce the overall uncertainty in the interferometer. 



B. radiation-pressure-fluctuations terms 



Next, we consider the contribution to {Az"^) from the radiation pressure fluctuations. 
This contribution is given by ( Ag^(t) )„: 



{Az')^,= {Aq\t))^ = {q')^-{q) 



2 

a 1 



(31) 



where the most dominant contribution is found to be 

pt pT p pt pr' p 

(Ag2)s = dr ds dAi\ / dr / ds' / dA\ 
Jo Jo J Jo Jo J 

' d^k r d^k' r d^~k r d^k' i i 



(27r)U (27r)U (27r)ij {27iyi^2k^2k' ^2k^2k' 
^^o^k(x|i, ^0, ■s)<9^o^k'(x||, 2;o, s) d^gUli^w, Zq, s')<9^o^ic'(*!l' ^0, -^0 



1 1 



+ ^2oWk'(x||, ^0, Wk(x||, 2;o, s) d^^Ul,m,zo, s')d^^U^^{±\\, zq, s') 



X 



Au A\\ 



lap / dr I ds I dr' I ds' 



'du 
2^ 



+ c.c. } + subdominant terms 

ZQ=L 



laPw A\\t^ Put 
' ' II cut = - 



Jo Jo Jo Jo 

+ subdominant terms 

3 



(32) 



m? (27r)3 

This result is valid for large time. All other subdominant contributions will be discussed 
in Appendix [b] The result of ( Ag^ at large time can be equivalently obtained by the 
particle-number operator approach, and it is given by the time integration of the velocity 
operator of the form 

(33) 



2w 

V = — n , 



m 
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where n is the particle number operator. In other words, the momentum received by the 
mirror can be effectively given by the momentum transfer of the scalar particles when they 
elastically bounce off the mirror. 

In contrast to the shot-noise type contribution that has the dependence, the stochas- 
tic motion of the mirror due to radiation pressure fluctuations yields an uncertainty pro- 
portional to the incident power P. This uncertainty can be attributed to the amplitude 
fluctuations of the incident field. These two types of noise have been extensively discussed 
in [H Hi [7] in the context of the interferometer. In particular, it is stated in [10] that shot 
noise given by the photon counting error can be understood in terms of phase fluctuations of 
the incident light while radiation pressure fluctuations are related to amplitude fluctuations. 
Our results are consistent with the arguments in [31 HI [3, (TU] . 

Other than contributions related to well-known shot noise and radiation-pressure fluctu- 
ations, modified field fluctuations also brings in a new component of noise. As the radiation 
pressure pushes the mirror into motion, the scalar field in the proximity of the mirror 
must change accordingly at every moment to enforce the boundary condition. In turn, this 
modification of field fluctuations will manifest themselves as quantum noise. Since these 
fluctuations arise naturally from the moving boundary condition of the fields, it must be 



included consistently in the expansion of Eq. (16), and this adds a new contribution into 
Az"^, as is described by (23). With the help of Eq. (19) and 



we have 

for ut ^ 1. Apparently, this term has mixed features that involve {q)a and {81!) a- Hence 
it may depend on both the incident power and the phase of incident state. In the end, 
however, the phase dependence cancels in this configuration so that the noise from modified 
field fluctuations depends only on the power of the input state. 

C. correlation terms 

In this framework, it is quite straightforward to compute the correlation between various 
noise. The correlation between shot noise and noise from random motion is given by 

( {A/ , Aq} )^^^ tan[u{L - zo)] , 

which can take either sign, depending on the distance between the detector and the mirror's 
original position. The presence of such sign-vary expressions gives the hope to reduce the 
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overall uncertainty ( A^;^ )„. It is worth to mention that the same result can be obtained by 



the particle- number operator approach using Eqs. (29), and (33). 

There exists another correlation term due to the additional source of noise, that is, the 
correlation between shot noise and the noise due to modified field fiuctuation, 

( { A/ , Add} ). - ^ tan[c.(L - zo)] , 

which has the same order of the magnitude as the the previous result. Furthermore, it leads 



to 



' ' [Add , Ag} + ( {A/ , Ag} + ( 1^/ , A/} 



Put^ /7 3 



— - tan [io{L — zq)] J , ^ 1 



m2 V2 2 

when we add up correlation terms of the order g^, as is summarized within the expression 



A'z),or inEq. (24). 



Now we may put all effects together, and the overall uncertainty {Az'^)^ is then given 

by 

(A.n.^^? + ^2C + ^fli-5c^V (35) 



Put 4 m ^ m2 \^ 2 2 ^ 
where ( = tan[c(}(L — zq)]. Interference between the incoming and reflected waves gives the 
quantity / the dependence on the effective distance between the mirror and detector. Later, 
the correlation between different sources of quantum noise can be found to possibly reduced 
the overall quantum noise. 

IV. EFFECT OF CORRELATION 

In the case of single-mode coherent state, C > represents positive correlation between 
shot noise and radiation-pressure fluctuations, while C < denotes negative correlation. The 
minimum value of ( Az'^ )a is found to be 



min(A.2), = C|2±^/^-^eU, (36) 



with an optimal value of the power P given by 



^°''* = =^v/22^^- ^^^^ 



In Eqs. (36) and (37), the "+ (-)" sign corresponds to the positive (negative) C case. 

The minimal value of {Az'^)a depends on (, but the range of the values of ( can not 
be arbitrary; instead, it has to be determined in a consistent way with the underlying 
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SQL 



FIG. 2: The minimal value of {Az'^)a as a function of ( in comparison with the SQL. The 
min ( Az^ is drawn in the unit of ^. 



assumptions. Slow motion of the mirror, -C 1, can be translated into an inequality, 



Poptwt^/m ^ 1, as can be seen from (19) and (28). The latter in turn implies \(\ < 1 by 
substituting the expression of Popt in (37) directly into the inequality. On the other hand. 



the definition of the effective displacement of the mirror in (17), formally by A/ ^ Az x dil, 
implicitly assumes 

\AzxdLl\>\l {Azfdll I . 
Together with the slow motion assumption, it gives 

2|CI 



uAz I < 



< 1 



The last inequality is imposed to ensure that the range of the value of Az is consistent 
with nonrelativistic motion of the mirror as well as the requirement on the definition of Az 
in (17), leading to \C,\ > \pl — 1 ~ 0.414. Thus |C| should lie within the range 



v^-l<|C|<l- 

The minimal value of ( Az^ as a function of C is shown in Fig. |2] 

For a positive min I^Az^ )a is a monotonic function of C, so it varies between 

1.78 — < min ( Az^ )s < 4 — , 
m m 



(38) 



compared with the SQL, Az'^^^^ = t/m. It is seen that we cannot beat SQL in the positive 



correlation case. For the upper (lower) bound in (38), the corresponding optimal power is 
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given by 

777- TTl 

Popt ~ 0.25 — and Popt ~ 0.09 — , 

respectively. On the other hand, for negative (, mm{A'^z)a is not a monotonic with (, 
having a maximum at ^ = —0.588, so in this case the largest value of min ( A^^;)^ is 

mm{Az'^)a ^ 0.136 — , 
m 

with the optimal input power 

0.132^. 

In the case of negative correlation, even the largest value of min ( A^z)^ can be lower than 
the standard quantum limit. Therefore, the value of the position uncertainty can be further 
reduced by other choices of (, and naively it can be as small as possible. However, there 



is a caveat. We have ignored higher-order terms in the series expansion (13) and in the 
calculation of ( Az"^ in terms of small q. In addition, we have employed the leading order 
approximation in the long-time limit by considering a large number of the particles in the 
coherent state. Hence min ( Az'^ )a cannot be arbitrarily small in that the subleading results 
will sooner or later contribute comparable effects. At any rate, we have shown that the 
standard quantum limit can be possibly overcome by establishing correlation between the 
intrinsic fluctuations of the field (shot noise) and its induced fluctuations due to the dynamics 
of the mirror (radiation pressure fluctuations and modified field fluctuations). 

So far, we have considered the single-mode coherent state. In reality, coherent state 
is at best prepared with a finite frequency bandwidth, so the observed response of the 
detector should be averaged over its bandwidth. In general, let the frequency distribution 
be described by some function of the form f{uj, u; do), where u is the mean frequency of the 
band and 2ao the bandwidth. The average over the distribution function is given by 

POO 

0{LD,ao)= du f{u,u]ao)0{u) . 



JO 

In our case, the quantity O is a product of a fast oscillating function of u due to the 
macroscopic scale L — zq and a relatively slowly varying component over the interval a) — ctq < 
u < u + cq. Suppose the bandwith is narrow in the sense that ctq ^ w and cro(L — zq) ^ 
1. Then the average can be approximated by substituting the variable u in the slowly 
varying component with the mean frequency u and directly performing average over the 
fast oscillating part. 

Thus if the coherent parameter a is independent of the frequency, the average of the 
normalization factor {8^1)^ is given by 

( dLl)l = A\a\^ sm^[2uj{L - zo)] x ~ 2\a\^ _1_ x , 
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from Eq. (15). Likewise, the corresponding frequency averages of the various sources of noise 



can be computed. The average of the shot noise term is given by 

and the average of the radiation-pressure fluctuations contribution is 

TAP),., = |«r siAuiL-zo)]cos^HL-zo)] x ^ \af x . 

It is interesting to note that the average of the rest of the terms totally cancels out by 
themselves, even though they contain drastically different information about correlation 
between noises and modified field fluctuations, 



\ ' rnf ' \ /cor 

The presence of finite frequency bandwidth reduces {AP{zQ,t; q + L)) to solely consist of 
contributions of shot noise and radiation-pressure fluctuations. Hence ( A2;2 )^ in the end is 
given by 

(A/2)^ _ 3 1 Poot^ 



(A^^)- = == = 7^ + ^' (39) 



which can be compared with Eq. (35). Minimization of {Az'^)a can be achieved with an 
optimal value of the power P, 

„ -\/3 m , ^ . 

Pop.-^^. (40) 



leading to the minimal value of {Az'^)a, 

min(A22\_ = v^l. (41) 
m 

The comparison will be made clear in Fig. [3] to highlight the significance of the correlation 
effects. 

At last we would like to make some remarks about the standard quantum limit in some 
other often-studied examples in the context of the measurement of the position of a quantum- 
mechanically free mass [2]. In those examples, the basic idea is that the evolution of the 
position of a free mass at time t is given by q{t) = q{0) + p{0)t/m. Thus, its position 
fluctuations with respect to a prescribed quantum state then turn out to be 



e t 



{Aq\t)) = {Aq\0)) + {Ap\0)) — + (Ag(O)Ap(O) + Ap(0)Ag(0)) - . (42) 

m 

It is assumed that the correlation term in the last expression is either vanishing or positive. 
The minimum position fluctuations is achieved by considering the zero correlation. Then 
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FIG. 3: a) Log-log plot of ( Az^ versus ^/Pt for various values of The straight line corresponds 
to the SQL. The parameters ^ = 10~^ and ^ = 10~^ are chosen, b) Log-log plot of ( Az^ 
versus y/Pt for various values of — . The straight line is the result of min( Az^)^. = 10""^ is 
chosen. 

by applying the minimum position-momentum uncertainty and minimizing the remaining 
expressions, we obtain the SQL. However, in [5], Yuen points out that the correlation term 
surely depends on the quantum state under consideration, and can be negative in principle. 
Thus, in general, beating the SQL is viable such that the position fluctuations at later time 
can be made small with a proper choice of the parameters, and even smaller than the SQL. 
The same conclusion that the SQL may not be a lowest limit of the measurement is shown 
in this mirror-field system. 

V. BACKREACTIONS INDUCED FROM RADIATION PRESSURE 

A stochastic description on dynamics of a mirror moving in quantum fields has been 
studied in [H]. There the associated semiclassical Langevin equation is derived with the 
method of Influence Functional by tracing out the quantum fields [H]. This Langevin 
equation incorporates not only the backreactions but also noises manifested by the field 
fluctuations. To examine the effect of the backreaction, let us write down the Langevin 
equation up to first order of the mirror's displacement, 

/■* SF 

mq{t)-^ J dt'&{t-t'){[Fo{t),Fo{t')])^q{t')-{Fo{t))^-{—{t))^q{t)=^^. (43) 

where 0(t) is the unit-step function, and the expectation value is taken with respect to the 
single-mode coherent state \a) of the scalar field. There are distinct contributions from the 
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quantum scalar field. In addition to the mean radiation pressure given by 

{Fo{t))^ = 11 dA^\ Af [z, t- z\ t] , (44) 



there are two different backreaction effects One is the local variation of the mean 

radiation pressure due to the mean displacement of the mirror, 

(^W).g(t) = \ jdA\\ dLAf[z,t-z',t]^^^,^^q{t). (45) 

The other comes from the nonlocal expression, associated with the retarded Green's function, 
which is constructed from the expectation value of the commutator of the forces. By the 
Wick's expansion, it can be written in terms of a product of the Wightman functions 
and xo of scalar field in the case of the single-mode coherent state, 

z e(t - t') ( [Fit), F{t')])^ = Af [z, t; z', t'] Xo [z, t; t'] ^^^^^^ , (46) 

with 

X^[z,t-z\t'] =ze{t-t') JdAiidA\ d,d,,{[<j){x),<j){x')])o. (47) 

Here we have used the point-splitting procedure to make this expression well-defined, and 
will set the small separation parameter e to zero in the end. In the coordinate representation, 
these two-point functions take the form 



Af{z,t; z',t') = {dMx)d,<Pix'))^ 

cos[cJ(2; — L)] cos[uj{z' — L)] | cos[cJ(t — t')] — cos[aJ(t + t' — 2L) — 2^?] | 



27r3 

1 



1 



t ( [0(x), (j){x')] )o = ^ [S{t -t'-R)- 6{t -t'-R)\-— \6{t -t' + R)- 5{t -t' + R) 
respectively. Here R and R are given by 

= (x|| - x'||)2 + {z- z'f , R^ = (x|| - x'||)2 + (z + z' - 2Lf , 

and in particular R describes the distance between the field point (xy , z) and the image of 
the source point (xy, z') with respect to the mirror at ^ = L. 

After plugging in the appropriate mode functions, the mean radiation pressure and its 
variation are given by 



F{t) )a = ^ ^^^11 { 1 - cos [2uj{t - L) - 2y?] } , (48) 
5q 



6F lap f ^ 

( —it) )^ = -y^^'^ll { sin[2a;(t - L) - 2^.] | . (49) 
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Although it is oscillatory, the mean radiation pressure of the incoming coherent state is 
always positive. Thus it keeps pushing the mirror along the positive z axis. Similarly we 



may carry out the integrals in Eq. (46), and end up with 



— I 



47r4 



co^\uj{z — L)] cos[cJ(2;' — L)] sin[cJ(t — L) — / dA\\dJl^\ 



(50) 



X 



b{t -t' - R) - ^S{t-t' - R) 
R 



J dt' e(t - t') sin [uj{t' -L)-ip] q{t') OA 
'A\\ sm\uj{t — L) — ipj |a;cos[a;(t — L) — 99] q{t) + sin[a7(t — L) — yj] g(t)| 



where we have taken the limit z 



L in the end. The outlines of the calculations can be 



found in Appendix C This nonlocal backreaction (50) reduces to local expressions. Espe- 



cially they include a term proportional to the first-order time derivative of the displacement. 
Although the phase of the incoming field changes with time, the coefficient of this q{t) term 
remains always positive. Therefore it plays the role of the frictional force, acting against 



the mirror's motion. It is seen that backreaction effects in Eqs. (49), (50) can be safely 



ignored as long as Put'^/m <^ 1 in the slow motion approximation, u}q{t) <C 1, as compared 



with the mean force term in Eq. (48). Some remarks are in order. These backreaction 



effects not only depend on the mirror's position and velocity but also have time-dependent 
coefficients. In particular, when the mirror is attached to an extra spring and undergoes 
oscillatory motion |27], the presence of a time-dependent coefficient of the position term may 
lead to instabilities on the mirror's dynamics due to parametric oscillation. Then, it might 
enhance the mirror's response of to small perturbations so as to improve the sensitivity 
for measurement of the weak signals [22] • Therefore the dynamics of the mirror, when the 
backreactions to the mirror are included, becomes rich. 

The accompanying noise to the backreaction from the commutator of the forces is denoted 



by in Eq- (43), which is associated with radiation pressure fluctuations. 



= Foit) - {Fo{t)), 



(51) 



as has been shown in Eq. (j5]). The consistent incorporation of these two effects in the 
Langevin equation relies on the existence of the relation between them, termed as the 
fluctuations-dissipation relation. For our mirror-fleld system in the case of single-mode 
coherent state, this relation can not be presented in a transparent way. Let us write the 
noise-noise autocorrelation function as 



1 

2* 



{m, m})^ = ^fb,t;z',t']ao[z,t;z',t'] 



z=L+e,z'=L ' 



(52) 
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where 



ao[z,t;z',t'] = ^JdAiidA\ d-A.' {{<p{x)A{x')])o. (53) 
is associated with anti-commutator of the scalar field in its vacuum state. This noise- 



noise autocorrelation function (52) can then be linked to the retarded Green's function 
in (46) through a relation between the Fourier transforms of the time-translationally invari- 
ant Green's functions xo and (Tq. When evaluated at the mirror surface z = L, z' = L, their 
Fourier transforms are respectively given by 

-0 = / S-o(-) (54) 

and the fiuctuations kernel (To{uj) can be shown to be related to the imaginary part of the 
Xo{^) kernel, 

ao{u) = Im Ixoioo)] [e{u) - e(-u;)] . (55) 

Thus we see that the role of the vacuum fiuctuations of the field seems to bridge the relation 
between the noise-noise autocorrelation function and the Green's function constructed by 
the commutator of the forces so that both effects can be treated in a consistent manner in 
the Langevin equation. 

The backreaction effect has been proposed theoretically in the work of Braginsky [23] 
for its possible role in setting the sensitivity limit of a detector in the Laser Interferometer 
Gravitational Wave Observatory (LIGO). Recently, radiation pressure cooling of a microme- 
chanical oscillator has been observed [21]. It raise a hope that quantum noise in a mirror-field 
system can be further reduced by properly tuning the backreaction effect. In Ref. [12], the 
backreaction has been extensively studied by modifying the dynamics of the test mass within 
the "optical spring" scheme to improve performance of the laser interferometer gravitational- 
wave detector. Here it is shown that these backreactions can be naturally incorporated in 
the Langevin equation by the field theoretic approach. Hence it will be a necessary next 



step that we apply the Langevin equation (43) to the interferometer in order to examine 
quantum noise as well as backreactions. 



VI. CONCLUSIONS 



The problem of quantum noise in the mirror-field system has been studied in the field 
theoretic approach. We consider that a single, perfectly refiecting mirror is illuminated by a 
single-mode coherent state of a massless scalar field, in addition to ambient vacuum fiuctu- 
ations. The net field is read out by a monopole detector, placed between the mirror and the 
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field source. The radiation pressure of tlie colierent state drives tlie mirror into motion. In 
tlie slow motion limit, we may identify different sources of quantum noise of the radiation 
field from readouts of the detector. In turn, the effective distance between the mirror and the 
detector can be obtained. The sources of quantum noise are respectively attributed to the 
intrinsic fluctuations of the field (shot noise), induced fluctuations arising from stochastic 
motion of the mirror due to the radiation pressure fluctuations, and modified field fluctua- 
tions which result from the mean displacement of the mirror. Their correlations can then be 
established resulting from interference between the incident field and the reflected field out 
of the mirror in the read-out measurement. The overall uncertainty can be found decreased 
(increased) due to negative (positive) correlation between noise sources. In particular, nega- 
tive correlation may lead to the situation that overall uncertainty is even lower than that in 
the standard quantum limit. Backreactions induced by the radiation pressure is studied by 
deriving the associated Langevin equation from first principles. The backreaction effects are 
found insignificant for a slowly moving mirror. The Langevin equation incorporates not only 
backreaction from radiation pressure on the mirror but also noise manifested by the field 
fluctuations. Equipped with the Langevin equation of the mirror-field system, it deserves 
further study to improve the performance of the interferometer as the work in [12] by taking 
advantage of the backreaction effect. 



Appendix A: Variation of / 



Since we have shown the to the order we are interested, the ordering of operators are 



irrelevant, we will just write Eq. (|13|) as 

/ 



In + qdLlo + ^q'dllo + 0{q^) 



LQ -t qULiQ -t -q u^iO -t ^{q ) , (Al) 

where the subscript denotes quantities evaluated on the mirror ai z = L. We further 



define AO = O — (O). Thus (Al ) up to first order in A becomes 



(/o) + A/o 


+ 


((g)+Ag)((9i/o)+A9i/o) 


1 r 
+ 2 


(/o) + A/o 


+ 


{q){dLlo) + Aq{dM + AdLlo{q) 



q) + Aq)\{dllo) + Adllo)'_ 



+ 



\qf{dlh) + 2Aq{q){dlh) + Adlhiqf + 0{A\ q') , (A2) 



Taking the expectation value yields 



{I) = {h) + {q){dM + ^{q)\dlh) + 



and thus the quantity A/ is given by 
A/ = A/o + 



1 



Aq{dM + AdLh{q)\ + [Aq{q){dlh) + -Adlh{q)' 



+ 



(A3) 
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This is Eq. (16). Since by definition of AO, we have ( (AO)^ ) = ( A^O ); hence the variation 
of / is 



A'I) = { A'lo ) + {A'q){dLlof + 2{Aq Addo ){q){dM + {A' {dd,) ){qf 



+ 



2{AhAq){dM + 2{AhAdM{q) 



(A4) 



+ 



2{AhAq){q){dlh) + 2{AhAdlh){qf 



where we have disregarded terms of the order q^ Aq and higher because the motion of the 
mirror is minute in comparison with other length scales like L, zq and o)"^. 



Appendix B: Evaluation of the subdominant terms in (A^g)^ 

The displacement q{t) of the mirror from its original position z = L can be expressed as 



q{t) = — I ds [ ds' I dA\\ : T^^ : (xn, z, s') 
mJo Jo J 



z=L 



(Bl) 



where the zz component of the normal-ordered stress tensor : T^z '■ (x||,-2,t) is given by 

: izz '■ [x\\,zA) = - lim OzO, / ^ , / ^ , 

X |Wk(x||, z,t)Wk'(xj|, z',t') afcOfc/ +Wk(x||,2;,t)Wk/(xj|,2;',t') 

+ Wk(X||, 2;,t)Wk'(X||) 4'^fc' +Wk(X||,Z,t)Wk/(xjj,z',t') 4/%| 

in terms of the mode functions Z//k(x||, z, t). It is convenient to have 

9,Wk(x||,z,t)|^^^ = (2^A;,)e^'^l|•'^l|-^-*e*'=^^ 
9,W^(x|,,z,t)|^^^ = (-2^A;,)e-^'^l|■^ll+*"*e-^^^^ 

and define 

/kk'(x!|,t) = 9,Wk(x||,z,t) X 9^.Wk'(x||,2;',t')|L~J=L 

= /k'k(^lh ^) 5 

5'kk'(x||,t) = dzUk{^\\,z,t) X (9^/Wk'(xj|, z',t')|2J^LL 

= _4yt k' ei(k||+k'|)-X||-i(i^+aj')tgi(fcz+fcOi 
= fi'k'k(X||,t) . 

The displacement operator q{t) then takes a simpler form 

— [ ds [ ds' (dA\\ f T- , / ; 

2m Jo Jo J {27c)lV2^J (27r)i 

The calculations of {A^q)^ involve computations of the expressions like 



/"* f f' d^k 1 f' d^k' 1 r t ^ 

k ^V'-^'i (2^71S/ Vv+h.c.| 
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k k' 



K«k'44')^ ^ («kak')«(44')"' 

('^kV44')^ - (4ak,)a(atat,)^, (aj^, 0^(3^4')" " («k'«k)'5(4'^L)°' ^^"^ complex 
conjugates, 



and so on. The result of the first term is given by Eq. (32). Here we would like to evaluate 



the second and the third terms, and compare them with the first term. 



1. Evaluation of second term 



Since 

«kak'44')" - («k«k')a(44)« = - - - + - ^)^(^' - 



+ 5{k - k)5(k - k)5{k' - k') + 5{k - k)5(k' - k)5(k' - 
+ 6{k - k)6{k' - k') + 6(k - k')6(k' - k) , 

and suppose that the particle number density |ap is sufficiently large such that the pure 
vacuum contribution can be ignored, the contribution from (akak'4'^L'^° ~ ('^k'^k')a(44''^'^ 
to the variance ( A^g(t) is then given by 



\a\ 
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|«P A 



ds j ds' jdA\\ dr dr' dA\^ 
Jo J Jo Jo J 

•^du u 



I ^3 



d^k I 1 1 
(27r)3 2^(27r)3 2^ 



^kk(x||,s')^kk(x||'^' 



m? (27r)3 2 Jo 27r 2 J^ Jq 







ds I ds' I dr I dr' e-^(^+^)^'e+*(^+^)"' 



(B2) 



where we have assumed that the incoming coherent field propagates along the z direction, 
i.e., k|| = 0. After integrating the time variables first, we are left with an cj-integral, which 
is formally UV-divergent. We regularize it with a cutoff frequency A, and then the result of 



(B2) is 



m^vr (27r)'^ 



In 



A + Cj 



1 + 



+ 



+ 0{ 



u \ A/ 

in the long time limit ut ^ 1. The cutoff-dependent term is logarithmic. For any sensible 
value of A, it should be of the order 0{1). Thus the expressions in the curly brackets are 



at most of the order unity. Compared with (32), this contribution is of about the order 
O{L0~^t~^) smaller; thus it is relatively subdominant. Note that the contribution from the 
expression {o\^ci\.iCirCi:r) a — {'Ai.'Ai') a is identically zero for a single-mode coherent 



field. 
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2. Evaluation of the third therm 



From 



k k' 



{aia^,alal)a - (4«k')a(44')'5 = '^(^ - k)5(k - k)5(k' - k') + 5(k - k)5(k' - k)5(k' - k) 



with a* being a complex conjugate of the complex coherent parameter a, its the contribution 
to the variance ( A^g(t) )„ is 



hi. ^^'/rf^ll/ dr[ dr'jdA', I ^^^^ /kk(x||,.')^^ak(xi|,r 



"^"^ '"ds Cds' tdr I dr' e-(-^)^'e+^(-+^)-' . 



'0 ^0 J JO Jo 

{271)3 2 2n 2j, j, j, j. 

We use the following approximations in the long time limit to simplify the cu-integral, 

fds fds' e-^-^y = 2 '^Ip^i^ _ 4_ ^ - -^^^ 

Jo Jo 



(B3) 



u — u 



(tU - W)2 



ZTT ; ^ d{UJ — 00) ~ I 



(CU — LU) 



t 



sm(u — uj)t 



Jo 



+ i 



u — u (w — uy 

t sin(a; + u!)t 



UJ + CO 



[00 + Uj) 



sinfa; + oo)t 



[u + oo {u + uy 

due to the fact that oot ^ 1 and oo > 0. Similar to the previous case, the integral over u is 



divergent, so a cutoff is introduced to regularize the result and Eq. (|B3|) becomes 

-2 



-i 2u}L 



ln^ + .^ + (9(i-; 
oo 4 oot 



m^vr (27r)^ 

There are three other similar terms and their results are either the same as this one or its 
complex conjugate. Following the arguments in the previous subsection, contributions from 



this type of expressions are negligible in comparison with (32) in the large time limit. We 



also note the contribution from (akak/a^ctj^,)a — ('2k'^k')'5('^^'^k/)a vanishes identically. 

As the final remarks, there are no contributions from (akOk/a^at,)^ — (a|^a|^,)a(a^at,)a 
and (ak'^k''^k'^k'^° ~ ('^k'^k')"('^k'^i{''^" because they vanishes identically, too. 



Appendix C: Evaluation of Eq. (50) 



We are interested in evaluating part of (50), denoted by J' 



J = J dAiidA\ J dt' e(t - t') sm[oo{t' - L) - ip] {d,d,> [(/>(x), (t){x)]^5q{t') 
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The retarded Green's function can be decomposed into two components; one corresponds 
to unbounded Minkowski space and the other to the image contribution with respect to 
mirror's original position z — L. We consider the unbounded component for the moment, 
and leave the image contribution to later sections. 

Let us take the partial differentiations dz before integrals over the cross sectional area 
^11 . Here we will assume the the linear dimension ^/A\\ is much larger than Uj~^. Thus the 
evaluation of may proceed as 

J= J dA\idA\\ dzdz'j dt' sm[iu{t' - L) - ip]9{t - t') S{t - t' - R) q{t') 

= J dA\\dA\ dzdz' |sin(i? - cui?) g(t - , -d = uj{t - L) - ^ , 



^ q^'^\t) / dA\\dAl\\ a,5,.{sin(i? - cui?) 7?""^} 

n=0 ' 



n=0 



°^ /_i\ri /-oo f 

2t^A\\ ^ {t) ^ dR jcu^e^i^'^-^ sin(i? - ujR) + u {R^ - e^) i?"-^ cos(i? - cui?) 

+ 2(n - l)a)e^i?"-^ cos(i? - GjR) - (n - 1) {R^ - e^) i?""^ sin(i? - 

- (n - 1) (n - 2)e2i?"-^ sin(i? - o^i?) i , 



where we have assume that the Taylor's expansion of q{t — R) with respect to R exists, and 
that implies the assumption that the memory of the system is rather short. We also have 
used the identities. 



dR _ _dR _z-z' _ e 

'dz ~ ~dz' ~ R ~ R' 
d^R d'^R 1 {z-z'Y 



R i?3 ' 



with e = z — z' ^ 0. For later convenience we will set e = |e|, and let £ — )■ in the end of 
the calculations. 

First, we check n = 2 case, and the integral over R becomes 



^ { 



ui'^e^ sin(i? - ujR) + a) ( - ^ ) cos(?? - uoR) 



2 / 2 

+ 2a) ^ cos(i? - uR) - ( 1 - ^ ) sin(i? - uR) 



uR cos{'& — uR) — sin(i? — uR) 
+ u'^e^sm{'d-QR) \ , 



+ 



uRcos{'d - uR) + sin(i? - a)i?)] 
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where 



dR 



[ 



i?2 

dR 



uR cos{d — uR) + sin(i? — cji?) j = e sin(i? — ue) , 
CuRcos{'d — uR) — sm(i? — = £sin(i? — ue) , 

/oo 
dR G?e^ sin(i? - UjR) = -a;^£ cos('j? - we) . 



Thus the integral /C2 is of the order 0[e) for the n = 2 case. Note that some of the integrals 
are still not well-defined in the large R regime so we introduce some cutoff to regularize 
them, but it turns out that the results are independent of the ultraviolet cutoff. 
Similarly for n = 1, we have 

/Ci = ^ dR jcj^^ sin(i? - ujR) + u - ^ cos{d - ujR) 



[ 



cue 



dR < ujcos('& — uR) + —r- u!Rsm('& — uR) — cos('& — uR) 

R^ 

dR ui cos{^ - uR) + 0{e) 



where 



dR 



cue"^ 



uR sin('i? — LoR) — cos('i? — uR) — —ue cos{'d — ue) . 



Finally in the case n = 0, the integral /Co is given by 



/Co 



dR 



u\^ 



^2 sin(i? -uR)+u\^- cos(i? - uR) - 2-^ sin(i? - uR) 
- 2'^ cos(^ - uR) +(^^-^\ sin(^ - uR) 



— —u cos(i? — ue) — —u cos ■& + 0{e) . 
In general, for n > 2, it can be shown that 



/Cn 



—ue cos{'d — ue) + n sin('i? — ue) 



so Kn behaves hke 0{e'^ ^) for n > 2. Therefore J is 

J = -2t:A\\ ^u cos d Sq{t) + sin ^ Sq{t) + 0(e) | 



(CI) 
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